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«N ; Abstract 

. We develop a method, based on Darboux' and Liouvihe's works, to find 

o , 

. first integrals and/or invariant manifolds for a physically relevant class of 

. dynamical systems, without making any assumption on these elements' form. 

We apply it to three dynamical systems: Lotka-Volterra, Lorenz and Rikitake. 



Q ; I. HISTORICAL OVERVIEW. 

^ [ In0i, Roger Liouville and A. Tresse developed a method for deciding whether two differ- 

ed I ential equations of the form 

d'^y f dy\'^ ( dy\'^ dy , . 

— + ai(x,i/) I — 1 +3a2{x,y)\—j + 3 as{x,y) — + ai{x,y) = (1) 

where the are arbitrary functions of the real or complex variables x and y, are geometrically 
equivalent, i.e. can be transformed into each other by the most general dependent and 
independent variable change 



x' = ip{x,y), y' = i){x,y) (2) 



PACS numbers: 02.30, 02.40, 02.90 
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This method was based on the construction of a "relative invariant" function called z/5 of the 
Oi and of their derivatives, such that in any transformation (|^) it becomes z/5 = J{x,y)~^ z/5 
where J{x, y) is the Jacobian of the transformation. In the general case, two equations such 
that their z/5 are non-zero and proportional to each other are indeed equivalent. If z/5 = for 
both, however, one cannot conclude at first, and other invariants, involving higher derivatives 
of the Oj, must be calculated in order to decide. As an application, Liouville proposed the 
effective reduction of Equation into its simplest canonical form, which in most cases 
leads to an explicit integration. 

Here we will adopt another point of view. We have derived from these theories a method 
for finding out first integrals for a wide and physically important class of dynamical systems 
without having to make any ansatz on their functional form. In the rest of this section, we 
shall recall some mathematical results of Darbouxi, Liouville and Tresse. Then we explain 
our method in Section and apply it in Section |lll| to three well-known dynamical systems. 



Finally, Section |^ discusses our results summarised in Table |I[ 



A. Essentials of Liouville theory. 

Consider a differential equation like (|l]). LiouvilleB defined the following functions — 

which are seen as functions of (x, y), forgetting the supposed relation between those variables. 

d ( dai \ d ( da^, da2 \ 

= Tr\^ 3 ai as + — — 2 — + ai 04 

ox \ox J oy \ oy ox j 

-3a2 ^^-2^ + aia4^ +ai ^^ + 80204^ (3) 

d f da4 \ d ( da^ da2 , \ 

-^1 = ^ + 3 02 04 - — 2 — — + ai 04 

oy \ oy J ox \ oy ox J 

( da3 da2 \ ( dai \ 

- 3 03 12 — - — + ai 04 I - 04 I - 3 ai ag I (4) 

/ dL2 dLA ( dLi dL2\ 

z/5 = L2 Li — L2 -T— ] + Li L2 ^ Li — — 

\ ox ox J \ oy oy J 

— ai Li + 3 a2 Li L2 — 3 a^ Li L2 + 04 L2 (5) 



The equation z/5 = means that 
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Li dx + L2 dy = (6) 

defines a particular solution of Equation (0). We shall call (^) a subequation of Equation (|lj), 
i.e. (|I|) is a differential consequence of (^). Notice that in the case L2 = the solution is not 
Li = — which would mean an unexpected lowering of this equation's differential order — 
but dx = 0, an absurdity as x is seen as the independent variable. Similarly if Li = 0, the 
solution is dy = 0, a solution possibly present in Equation (|l|) but not very interesting. 

Now suppose neither Li nor L2 vanish identically and define a = —L2/L1. Then 

f da da s o 2 o \ /«n 

i/5 = a — + — + ai a + 3 02 a + 3 as a + 04 1 [7) 

and Equation (||) can be rewritten dy/dx = a{x,y). Conversely suppose there is a first- 
order subequation to Equation (|1|), namely dy/dx = A{x,y). Then A is a solution to the 
first-order non-linear PDE 

dA dA 

A — + — + aiA^ + 3a2A^ + 3a3A + a4 = (8) 
oy ox 

we shall discuss later. Liouville theory is a finite effort tool for finding particular solutions 
to Equation (||). 

B. Darboux polynomials and first integrals for polynomial dynamical systems. 

Consider an autonomous polynomial dynamical system 

Xi = Vi{x), i = l...n (9) 

We say that a polynomial /(xi, . . . is a Darboux polynomial of if there exists a 
polynomial "eigenvalue" p such that 

dt {r[ * dxi 

In other words, there is an algebraic variety, defined by /(xi, . . . , Xn) = 0, which is invariant 
by the flow of V . In this respect, this notion is a neighbour of the notion of subequation we 
have seen. 
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Darboux polynomials are tools for finding outiH, but also proving the non-existence (cf.i 
for an example) of first integrals to polynomial dynamical systems. We shall not enter into 
the details. Let us just notice that a polynomial first integral is simply a Darboux polynomial 
with eigenvalue 0; and that a Darboux polynomial / with constant eigenvalue a gives rise 
to the time- dependent first integral /e~"*. More rational and algebraic first integrals can 
be built with the "basic blocks" of Darboux polynomialsi'i; and, conversely, a theorem of 
Brun^ says there cannot be an algebraic first integral of (P) unless there is a rational one, 
which in turn implies the existence of Darboux polynomials. 

In brief, the problems of existence of first integrals and Darboux polynomials for a 
polynomial dynamical system are very tightly linked. Notice also that all these objects, like 
ordinary eigenvectors and -values of linear endomorphisms, naturally live in C. 



II. PRINCIPLE OF THE METHOD. 

We shall draw our interest to autonomous three-dimensional polynomial dynamical sys- 
tems which are of first degree in one of their three variables, e.g. z. Their general form is 
thus: 

'x= Vl{x,y) + zVl{x,y) 

< y= V^y{x,y) + zVl{x,y) (10) 
, i = V%x,y) + zVl{x,y) 

which we may abbreviate as X = V{X). Dynamical systems of this kind are frequently met 
in physics: well-known examples are the Lorenz model, or the various three-wave interaction 
problems (Rabinovich etc.). Very often they are indeed of first degree in all their variables. 
We can use this feature for harvesting more information — an example is given in the 
paragraph about the Lorenz model. 

We assume to have found out and studied all solutions with a; = est. Assuming x noncon- 
stant, we shall transform the system (|T^) into a non- autonomous second-order differential 
equation linking y and x which will turn out to be of type (|l]). 



FIRST INTEGRALS FROM GEOMETRICAL EQUIVALENGE 5 

Now we settle in a region of space where x ^ and take x as the independent variable, 
parametrising the integral curves of (0). The relation 

{Vlix, y) + z Vlix, y)) ^ = V^x, y) + z V\{x, y) (11) 

is satisfied along all integral curves. Hence, writing p = dy/dx, 

z{vlp-Vl)=Vl-pVl (12) 

These equations define the mappings (f):{x,y,z) i— {x,y,p) and (j)*: {x,y,p) ^— {x,y,z), 
which are homographic and hence: 

1. They are one-to-one wherever they are defined and their determinant C{x,y) = 

~ is non-zero; as it involves only the variables x and y the surface 

S = {C{x,y) = 0} can be seen either as a submanifold in the {x,y,z) space or in 
the (x, y,p) space. 

2. The surfaces 5*1 = {V^{x,y) + zVl.{x,y) = 0} in the {x,y,z) space and 82 = 
^Vl.{x,y) p — Vy{x,y) = o| in the {x,y,p) space are singular. Any point on S*! \ S 
is sent to p = 00: this happens when x = 0, and the tangent to the integral curve is 
orthogonal to the x-axis, i.e. "vertical" in {x, y) representation {dy/dx = 00). Similarly 
any point on 6*2 \ S is sent to 2; = cxd. 

3. On S, and 0* are "constant along fibres", i.e. two points of S having different z 
(or p) are sent to the same image, having the same {x,y) as the original point, hence 
lying on S. Thus it is always possible to get it as the image of a point on S \ 6*1 and 
calculating it that way shows that 0(S) = S fl 5*2; and similarly 0*(S) = S fl ^i. 

Since we are concerned with a differential problem, we have to study what the vector field 
V{X) becomes under the action of the tangent map Tx4>- And, indeed, points X on S 
differing only in the z coordinate have the same image by (p, but different V{X) such that 
the corresponding Tx(f){V{X)) also generally differ. As all these vectors are attached to the 
common image of the points X, this can cause a loss of information, leading, as we shall 
see, to important practical difficulties. 
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Similarly, we find that 

{vUx, y) + z Vl{x, y)) ^ = Vl{x, y) + z V\{x, y) 

We calculate dz/dx by differentiating (|12|) with respect to x, putting the result into the 
previous formula, and then replacing z itself with its value in function of p given by Equa- 
tion (|r^). This leads to a differential equation in p which, in Cauchy form, reads 

dp _ N{x,y,p) 

dx C{x,yy ^ ' 

where is polynomial in {x,y,p) and of degree three in p. Interpreting p as dy/dx, we see 
Equation (|1^) as a differential equation of Liouville type like (|I|). There are two essential 
facts in these computations. One is that the denominator of ([13D is exactly the square of the 
determinant C{x,y) of (p, so Equation ([T3| ) will not set any further problem as long as its 
construction is valid. The other one is that Equations ([T0|) are of degree one in z: it ensures 
not only the good behaviour of the z ^ p correspondence but also the Liouville form of the 
differential equation (|13D . 

We intend to apply Liouville theory to Equation (|13]) in order to obtain subequations 
for it. Now, we must take care of their possible relationships with the forbidden surfaces. If 
a subequation defines a curve in the {x,y,p) space which is not contained in S or 5*2, there 
is no problem: it will be pulled back into the (x, y, z) space by the 0* map, which coincides 
then with the reciprocal of (p. 

But in the computation of Equation (p!3D, we have used the z = 4>*{x, y,p) map, and then 
suppressed the denominator V\p — V\^. Thus, the singular manifold at 5*2 has disappeared 
in (|I^). But consider a curve plotted on 5*2 (i.e. V\.p — V]^ = which is, moreover, a jet 
(i.e. p = dy/dx). Then, identically 

dy ^ ^ Vl{x,y) 

dx ' Vl{x,y) 

We can check that this A is always a solution to Equation (^), whatever the vector field V 
may be. Hence any jet plotted on S2 is a subequation of Equation j^T^j). However, this jet 
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cannot yield an invariant manifold in the (x, y, z) space unless it is made of images by of 
points in this space. Now, it is easy to show that the only points on 5*2 that can be written 
as (f){x,y, z) are those on S fl 5*2; therefore, as we have seen, they are images of points also 
lying on S. 

Yet we know that the puUback of the vector field V is not necessarily well-behaved on 
S. Thus, the system (0) and the differential equation (13) can behave quite differently on 
S. 

If we find as Equation (|^) the equation of a jet on S'2, we have to check independently 
whether S is an invariant manifold for V or not. In dynamical systems containing parame- 
ters, this can be rephrased as: find at what condition on the system 's parameters the equation 
denominator C{x, y) is a Darhoux polynomial for the system PT[/ 

III. RESULTS FOR SEVERAL DYNAMICAL SYSTEMS. 

We have applied the method exposed in Section |l| to three different dynamical systems of 
type ([T0|) depending on real parameters: Lotka-Volterra, Lorenz, Rikitake. We shall discuss 
the results obtained in the rest of this section. 

A. The (a, 6, c) Lotka Volterra system. 

This remarkably symmetric system 

' X = X {cy + z) 

y = y{az + x) (14) 

z = z{hx + y) 

appeared first as a model for three-species competition, yet has been found later in plasma 
physics. A considerable amount of research has been done on it, using many techniquesHiHl^. 
Here we shall follow the process exposed in Section |I[ 

Since Equations (|Iip are invariant by simultaneous circular permutations of (x, y, z) and 
(a, 6, c), it is equivalent to perform the method with any couple of variables. Once this is 



FIRST INTEGRALS FROM GEOMETRICAL EQUIVALENGE 8 

done, more results can be got by the above symmetry. There is also a symmetry in taking 
X = x'/b, y = z'/c, z = y'/a and a = 1/c', b = 1/6', c = 1/a', which will appear in the 
distribution of the z/5 = cases. 

So, we take z as the independent variable and eliminate x, and find 

^ ^ {b-l){l + abc)Ql,{y,z) 
^ z{y~abzY ^ ^ 

and 

_ {b-l){l + abc)Ql,^{y,z) 

z{y-abzy ^^^^ 

and 

_ ib-lf{l + abcfP,Uy,^) 

z^y^{y-abzY^ ^ ' 

where Pabc^ Qlbc^ Qabc polynomials whose coefficients depend on (a, b, c) and which we 
do not write down for the sake of brevity. The cases 1 + abc = and 6=1 are known: 
the first one is the full integrability of the system, with, in particular, the first integral 
abx + y — az; in the second one we have the Darboux polynomial y — az whose eigenvalue 
is X. We remark this Darboux polynomial is exactly the denominator. . . Notice also that in 
those two cases, Li and L2 vanish together with z/5 so that Equation @ is an identity and 
cannot be used for finding out Darboux polynomials. 

Now the cases where all coefficients of Pabc are zero are listed below in Table |I[ We 
notice the presence of the symmetry a = 1/c', b = 1/6', c = 1/a' in this list; one case 
(a = 1, 6=1, c = 1) is self-symmetric. We shall handle in some detail one of these "exotic" 
cases, viz. (a = 1/4, 6 = 2, c = —5). The subequation (D reads: 

y (^-3 z'^ + 16 zy- 32 y^) (z^ - 7 zy + IQy"^) =0 (18) 

Reverting to the original variables (x, z) changes this equation in 

y{z-2y){lQy'^ -2xz-^yz + z^)=Q (19) 
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The expression z — 2y is proportional to the equation denominator; its presence here is 
an artefact due, as we have seen, to a former suppression of denominator. It should not 
be taken in account since it corresponds to b = 1. On the other hand, the other two 
factors are Darboux polynomials, since their derivatives with respect to the system (|H|) are 
y = y {z/A + x) and 

(iQy'^ - 2x z - 8y z + z^) = 2x (iQy^ - 2x z - 8y z + z^) 

This illustrates the validity of the method in the general case. The results for all cases 
are summarised in Table |. We get no information for the self-symmetric case since it is a 
specialisation of 6 = 1, so our method cannot be applied. 

B. The Lorenz model. 

Another well-known and intensively studied&ii dynamical system, the Lorenz model 

X = a iy — x) 

y = Qx — y — xz (20) 

z = xy — b z 

originally thought as a simple model for atmospheric turbulence, was the first example 
of a low-dimensional chaotic deterministic dynamical syste m0. All known first integrals 
have been obtained or reobtained by Kusi, using the non-decisive procedure of Carleman 
embedding. Here we shall recover some of them methodically by the means of Liouville 
theory. 

Since there is no symmetry among variables here, we can proceed three times to the 
calculation of the L12 and 1/5 functions, eliminating each time one of the three variables. 
Indeed, the Liouville equation like (|l^) contains no more information than the dynamical 
system like (0) does, but it has more singularities; yet, as we have seen, these singularities 
often contain useful matter about the dynamical system's invariant manifold structure. 
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Eliminating z and choosing x as the independent variable yields C{x,y) = a x {y — x). 
Neither x = nor y — x = can be interesting invariant manifolds, since both imply x = 0, 
so X = est and y = est. This, in turn, also implies that z is constant, and the manifold 
reduces to a fixed point. 

The functions L2 and z/5 are (cf.0) 

l + b + a 



a{y-xf 



(21) 



and 



= {l + h + a)P,^,{x^y) 

In the obvious case 1 + h + a = 0, we also have L2 = 0. Thus, as explained 
in Section |, our choice of variables was a bad one. Looking for other cases, we 
only get three points in the (b,a,g) space, viz. (6 = 0, cr = — 1), (6 = 2/3,(7 = 1/3) 
and {b = —16/5, a = —1/5, g = —7/5). They are specialisations either of known caseslli 
or of the 1 + b + a = case. In those cases, we get subequations that do not give rise to 
Darboux polynomials, i.e. equations that represent jets plotted on the surface 82. 

Now, if we eliminate y and choose x as the independent variable, we get as equation 
denominator C{x, z) = b z — x"^. Its derivative with respect to the system (^0]) is 

^(bz-x^) ={b-2a)xy-(b^ z-2ax^) (23) 
= -2(x{bz-x'^) + {b-2cx) (xy-bz) 

The remainder is of degree one in x. Thus C is a Darboux polynomial iff 6 = 2 cr; then the 
eigenvalue is —2 a, so / = (a;^ — 2 a z) e^°"* is a first integraJii. 

We have calculated 1/5 and found 

_ ib-2a){l + b + a)Pf,., 

a^bz-x^y^ ^^^^ 

PbcTQ being such that its coefficients never simultaneously vanish. In the case b = 2a, 
Equation (|^) yields x — a dz/dx = 0. This is the equation of a jet on 82, and since we are 
in the good case, we find S as the invariant surface. 
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When 1 + 6 + 0" = 0, Equation still yields x — a dz/dx = 0. But in this case, S is not 
invariant and we do not have a Darboux polynomial. 

Finally, we have taken z as the independent variable and eliminated x. We find C{y, z) = 
b gz — bz"^ — y"^ and 
dC 

— = xy {{b-2) Q + 2{l-b)z) + 2y'^ + 2b^ z^ - b^ Qz (25) 
= -2{bQz -bz"^ - y'^) + R{x,y,z) 

the remainder R being of first degree in y; hence C is a Darboux polynomial iff 6 = 1 and 
g = 0. Then, dC/dt = —2C and we get that way the first integral / = (z^ + y'^) e^*0. 



As for z/5, it is equal to 



(l + b + a)yPb„g{y, z) 
[bgz — bz^ — y^)^^ 



where Pbag{y,z) = iff 6 = 1 and ^ = 0. Let us handle first the latter case. In that 
case, Li = —azy/^z"^ + y'^Y and L2 = —ay'^/{z^ + y^)^, so Equation (|^) simplifies as 
z dz + y dy = 0. This is the equation of a jet on S2, and we get S as invariant manifold. 

Now, in the case 1 + b + a = 0, there is another simplification in Equation (^), namely 
y Li = {z — g) L2 and hence {z — g) dz + y dy = 0. But this is the jet on 5*2, so we get no 
information in this case. 



C. The Rikitake dynamo. 

This dynamical systemi 

'x = —fix + y {z + j3) 

< y = -fiy + x{z-p) (27) 
i = a — xy 
models the variation of the earth's magnetic field with time. 

Let us take x as the privileged variable and eliminate z. The denominator in Equa- 
tion (0) is C{x, y) = fi {y'^ — x"^) + 2 (3 X y, and 
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where the coefficients of Pa/s^i cannot vanish together unless /3 = or /i = 0. 
The derivative of C{x,y) with respect to the system (|27|) is 

-^ = (y^-x'){2(3'-2fi')-2/3xy + 2Pz{x + y) (29) 

Assume C is a Darboux polynomial of eigenvalue P{x, y, z) = A{x, y) + z B{x, y). Then the 
identification of the z terms in Equation (^) gives B C = 2[3 {x + y). Since C is of second 
degree, this is impossible unless /5 = 0. 

When /3 = 0, Equation reads dC/dt = — 2/iC, so C is a Darboux polynomial of 
this system, which gives the first integral (y^ — x^) e^'^*. On the other hand, 

hence Equation becomes —xdx + ydy = 0. This is a jet on 5*2, but we are in the "good" 
case, and we recover the Darboux polynomial y"^ — x"^. 

If now /i = then 

_ -a{3x^ + 2x^y^ + 3y^) _ a {3 x^ + 2 x^ y^ + 3 y^) 

-^1 - . no A ) -^2 



4^23.3^4 ' AP^x^y^ 

so the subequation is once more —xdx + y dy = 0, the jet on 82- Hence we do not obtain 
any Darboux polynomial unless f3 = 0. 

We have also performed the computations with the other two couples of variables. They 
have not given more information than the previous ones. 



IV. CONCLUSION. 



We have obtained, by a methodic procedure, numerous cases of Darboux polynomials for 
the Lotka-Volterra system. This "Darboux polynomial searcher" can be seen as an input 
to algorithms which need Darboux polynomials, such as the Prelle-Singer procedure. Up 
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to now, that procedure began with a systematic search, which obhged to set an a priori 
hmit on the polynomial's degree in all its variable^l. Some new results! allow to refine the 
search by restricting the choice of the possible highest-degree homogeneous components of 
the tentative Darboux polynomials, while speeding it up when the system's coefficients are 
rational numbers. Yet they are valid for dynamical systems of dimension 2, and until now 
have no counterpart in dimension 3. 

Our method has also reobtained some known first integrals for the Lorenz and Rikitake 
systems, though all cases have not been found, and despite the "divergence enigma" we now 
explain. 

In both Lorenz and Rikitake systems, the divergence is a constant, respectively —1 — b — a 
and —2 /x. In both cases, its vanishing triggers the vanishing of z/5, but also the reduction 
of Equation (^) to a singularity from which no information can be extracted. However, the 
Rikitake dynamo possesses when /x = a time-dependent first integral I = x'^ — y'^ + 4:f3z — 
Aa f3t; this kind of first integral cannot be detected by our method, since it does not arise 
from a Darboux polynomial, but from a polynomial f{x,y,z) such that df/dt = est. Such 
a first integral may exist only for dynamical systems having a constant term, so there is no 
chance to find any for the Lorenz model. But there may be a first integral of some special 
kind — indeed, numerical experiments exhibit a regular behaviour when 1 + 6 + a = 0. 

Table | summarises the results obtained by our method. Abbreviations are: DS for 
"dynamical system", FI for "first integral", and DE for "Darboux element", i.e. a couple of 
polynomials (/, p) such that df /dt = p f. 



FIRST INTEGRALS FROM GEOMETRICAL EQUIVALENCE 



14 



REFERENCES 

^ R. Liouville, Memoire sur les invariants de certaines equations differentielles et sur leurs 
applications, Journ. Ec. Poly. 49 (1889), 7-76. 

^ A. Tresse, Sur les invariants differentiels des groupes continus de transformations, Acta 
math. 18 (1894) 1-88. 

^ G. Darboux, Memoire sur les equations differentielles algebriques du premier ordre et du 
premier degre, Bull. Sc. Math. 2eme serie, 2 (1878) 60-96, 123-144, 151-200. 

^ Y.K. Man, First integrals of autonomous systems of differential equations and the Prelle- 
Singer procedure, J. Phys. A 27 (1994) L327-L332. 

^ Y.K. Man, M.A.H. MacCallum, A Rational Approach to the Prelle-Singer Algorithm, to 
appear in J. Symbolic Computation. 

^ S. Labrunie, On the polynomial Grst integrals of the (a, b, c) Lotka-Volterra system, to 
appear in J. Math. Phys. 

^M. Kummer, R.C. Churchill, D.L. Rod, On a result of Bruns, Can. Math. Bull. 33 (1990) 
175-180. 

^T.C. Bountis, A. Ramani, B. Grammaticos, B. Dorizzi, On the complete and partial 
integrability of non-Hamiltonian systems, Physica A 128 (1984) 268-288. 

^ J.M. Strelcyn, S. Wojciechowski, A method of finding integrals for three-dimensional dy- 
namical systems, Phys. Lett. A 133 (1988) 207-212. 

°B. Grammaticos, J. Mouhn-OUagnier, A. Ramani, J.M. Strelcyn, S. Wojciechowski, In- 
tegrals of quadratic ordinary differential equations in R^: the Lotka-Volterra system, 
Physica A 163 (1990) 683-722. 

^E.N. Lorenz, Deterministic Nonperiodic Flow, J. Atm. Sci. 20 (1963), 130-141. 



FIRST INTEGRALS FROM GEOMETRICAL EQUIVALENCE 15 

1982 H. Segur, Solitons and the Inverse Scattering Transform, in A.R. Osborne and 
P. Malanotte Rizzoli (eds), Topics in Ocean Physics, North- Holland, 1982, pp. 235-277. 

^^M. Kus, Integrals of motion for the Lorenz system, J. Phys. A 16 (1983), L689-L691. 

G. Levine, M. Tabor, Integrating the nonintegrable: analytic structure of the Lorenz 
system revisited, Physica D 33 (1988), 189-210. 

^^T. Sen, M. Tabor, Lie symmetries of the Lorenz model, Physica D 44 (1990), 313-339. 

V.S. Dryuma, Projective duality in the theory of second-order differential equations (in 
Russian), Mat. Issled. (Kishinev) 112 (1990) 93-103. 



FIRST INTEGRALS FROM GEOMETRICAL EQUIVALENGE 



16 



TABLES 



DS 




Denominator 


Parameters 


Information obtained 








(0,2,c) 


DE: {y;x) 








(1/(2 6), 6,1) 


DE: {bx- z;y) 








(1,6,2/6) 


DE: {x — cy; z) 








(1/4,2,-5) 


DE: (16y2-2xz-8yz + z2;2x) 


Lotka 


(y,^):(0) 


y — abz 


(-1/5,1/2,4) 


DE: (100y2 - 2bxy + ^y z + 4.z'^;x) 


{a,b,c) 






(1,2,-2) 


DE: (y2 — x z — y z;2x + z) 








(-1/2,1/2,1) 


DE: (-2 xy + 2y2; + z^;x + y) 








(1,1,1) 


none 








(-1/2,0,1) 


none 








(0,2/3,1) 


DE: (3z - 2x;y) 




other 


cycUc permutation of the above results 








(-1 - a, a, g) 


none 






x-y 


{2/3,1/3, g) 


none 








(-16/5,-1/5,-7/5) 


none 


Lorenz 




hz — 


(2 a, a, g) 


FI: (x2 - 2f7 2)e2<^* 


(6, cr, g) 






(-1 - cr, a, g) 


none 






bz{g- z) -y^ 


(l,a,0) 


FI: (z2+y2)e2t 








(-1 - a, a, g) 


none 


Rikitake 




H (y^ - x^) 


(a,0,/x) 


FI: (y2 _^2)g2Mt 


(a,/3,/x) 




+2Pxy 


(a,/3,0) 


none 




other 




nothing more 



TABLE I. Results obtained by our method 



